Enhance synchronizability by structural perturbations 
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In this paper, we investigate the collective synchronization of system of coupled oscillators on 
Barabasi- Albert scale-free network. We propose an approach of structural perturbations aiming at 
those nodes with maximal betweenness. This method can markedly enhance the network synchro- 
nizability, and is easy to be realized. The simulation results show that the eigenratio will sharply 
decrease to its half when only 0.6% of those hub nodes are under 3-division processes when network 
size N = 2000. In addition, the present study also provides a theoretical evidence that the maximal 
betweenness plays a main role in network synchronization. 
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I. INTRODUCTION 



Many social, biological and communication systems 
can be properly described as complex networks with 
nodes representing individuals or organizations and edges 
mimicking the interactions among themHHg. One of 
the ultimate goals of the current studies on topological 
structure of networks is to understand and explain the 
workings of systems built upon those networks: for in- 
stance, to understand how the topology of Internet affects 
the spread of computer virus |J, |fj, |fj, how the structure 
of power grids affects the cascading behavior0,@,@; now 
the connecting pattern of intercommunication networks 
affects traffic dynamics [Tol ITU fl2| . and so on. 

In the past few years, with the computerization of data 
acquisition process and the availability of high computing 
powers, scientists have found that many real-life networks 
share some common statistic characteristics, the most im- 
portant of which are called small- world effect milil 

and 

scale- free property 0]. The recognition of small- world 
effect involves two facts, small average distance fas 
L ~ InJV, where N is the number of nodes in the net- 
work) and large clustering coefficient 01 • The number of 
edges incident from a node x is called the degree of x. 
The scale-free property means the degree distribution of 
the network obeys power-law form, that is p(k) ~ k 1 , 
where k is the degree and p(k) is the probability den- 
sity function for the degree distribution. 7 is called the 
power-law exponent, and usually between 2 and 3 in real 
world[Hll0]. This power-law distribution falls off much 
more gradually than an exponential one, allowing for a 
few nodes of very large degree to exist. Networks with 
power-law degree distribution are referred to as scale-free 
networks, although one can and usually does have scales 
present in other network properties. 

Synchronization is observed in a variety of natural, so- 
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cial, physical and biological systems [Tsl Il9l l2pt l2lj| . To 

understand how the network structure affects the syn- 
chronizability of the network is of not only theoretical 
interest, but also practical value. There are many previ- 
ous studies about collective synchronization, with a ba- 
sic assumption that the dynamic system of coupled os- 
cillators evolves either on regular networks [2^ |23L |24| . 
or on random ones[2^, [2(|. Very recently, scientist fo- 
cus on synchronization on complex networks, and find 
that the networks of small-world effect and scale-free 
property rnay be easier to synchronize than regular 
lattices |H El El US HI IH . 

Since there are countless topological characters for net- 
works, a natural question is addressed: what is the most 
important factor by which the synchroizability of the net- 
works is mainly determined? some previous works indi- 
cated the average distance L is one of the key factors, the 
smaller L will lead to better synchronizability psl I2H l3l| . 
Other researchers focus on the role degree heterogene- 
ity plays. They found greater heterogeneity will re- 
sult in poorer synchronizability, and demonstrated that 
the maximal betweenness |33t l34j B max may be a proper 
quantity to estimate the network synchronizability. With 
smaller Bmnj 1 the network synchronizability will be 
better|3!| [36|. However, the above results and conclu- 
sions are still debated. 

In this paper, we investigate the collective synchroniza- 
tion of system of coupled oscillators on Barabasi- Albert 
scale- free networks (B A networks for short) We pro- 
pose an approach of structural perturbations, which can 
markedly enhance the network synchronizability, and is 
easy to be realized. It also provides a theoretical evi- 
dence that the maximal betweenness plays a main role in 
network synchronization. 

This paper is organized as follow: in section 2, the 
conception of synchronizability will be briefly introduced. 
In section 3, we will describe the approach of structural 
perturbations. Next, the simulation results will be given. 
Finally, in section 5, the conclusion is drawn, and the 
relevance of this approach to some real-life problems is 
discussed. 
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FIG. 1: The sketch maps for 3-division process on xo- The 
solid circle in left side is the node xo with degree 6. After 
3-division process, the xo is divided into 3 nodes xo,xi and 
X2 that are fully connected. The six edges incident from xo 
redistribute over these three nodes. 

II. SYNCHRONIZ ABILITY 

We introduce a generic model of coupled oscillators on 
networks and the master stability function|37j|. which is 
often used to test the stability of the complete synchro- 
nized stats. 

Each node of a network is located an oscillator, a link 
connecting two nodes represents coupling between the 
two oscillators. The state of the ith oscillator is described 
by x l , we get the set of equations of motion governing the 
dynamics of the N coupled oscillators: 

JV 

x'^Ffx'J+a^W), (1) 
j=i 

where x l = F(x*) governs the dynamics of individual os- 
cillator, H(x 3 ) is the output function and a the coupling 
strength. The JV x N coupling matrix G is given by 

!— fcj for i = j 
1 for j 6 A,. (2) 
otherwise 

All the eigenvalues of matrix G are nonpositive reals 
because G are negative semidefinite, and the biggest 
eigenvalue 70 is always zero because the rows of G 
have zero sum. Thus, the eigenvalues can be ranked as 
7o > 7i > ' ' ' > TiV-ij and the synchronization manifold 
is an invariant manifold, that is, the fully synchronized 
state: x 1 = x 2 = • • • = x N = s satisfies s = P(s). 

Let £ l be the variation on the ith node, and the col- 
lection of variation is £ = (t; 1 ,^ 2 , ■ ■ ■ = £ N ), we get the 
variational equation of (1) 

e= [ljv O DP + o-G ® DH]f. (3) 

Diagonalizing G in the second term of (3), a block diag- 
onalized variational equation is obtained and each block 
has the form 

ik = [DF + tnkD^ k , (4) 

where D denotes the Jacobian matrix, 7^ is an eigenvalue 
of G, k = 0, 1, 2, • • • , JV — 1. k — corresponds to the 



mode that is parallel to the synchronization manifold. 
Let a = <77fc, rewrite (4) as 

C= [DF + aDH}{. (5) 

Since DF and DH are the same for each block, the largest 
Lyapunov exponent A max of Equ. (5) only depends on a. 
The function A max (a) is named master stability function, 
whose sign indicates the stability of the mode: the syn- 
chronized state is stable if A max (o!) < for all blocks. 

For many dynamical systems, the master stability 
function is negative in a single, finite interval ((11,02) 
the largest Lyapunov exponent is negtive|33. Therefore, 
the network is synchronizable for some a when the eigen- 
ratio r = 7jv-i/7i satisfies 

r < a 2 /ai (6) 

The right-hand side of this equation depends only on the 
dynamics of individual oscillator and the output func- 
tion, while the eigenratio r depends only on the coupling 
matrix G. The problem of synchronization is then di- 
vided into two parts: choosing suitable parameters of 
dynamics to broaden the interval (ai, 0.2) and the analy- 
sis of eigenratio of the coupling matrix. The eigenratio r 
of coupling matrix indicates the synchronizability of the 
network, the smaller it is the better the synchronizability 
and vice versa. In this paper, for universality, we will not 
address a particular dynamical system, but concentrate 
on how the network topology affects eigenratio r. 

III. STRUCTURAL PERTURBATIONS 

As mentioned above, the nodes with very large be- 
tweenness, namly hubs, may reduce the network synchro- 
nizability. So the present method of structural perturba- 
tions aims at these hubs. For a hub xq, we add m — 1 
assistant nodes around it, labelled by x±,X2, ■ ■ ■ ,x m —\. 
These m nodes are fully connected. Then, all the edges 
incident from xq will relink to a random picked node 
Xi(i — 0, 1, ■ • • , m — 1). After this process, the between- 
ness of xq is divided into m almost equal parts associating 
with those m nodes. We call this process m-division for 
short. A sketch map of a 3-division process on node £0 
is shown in Fig. 1. 

Due to the huge size of many real-life networks, it 
is usually impossible to obtain the nodes' betweenness. 
Fortunately, previous studies showed that there exists 
strongly positive correlation between degree and be- 
tweenness in BA networks and some other real hetero- 
geneity networksjl^, that is to say, the node with 
larger degree will statistically have higher betweenness. 
Therefore, for practical reason, we assume the node with 
higher betweenness is surely of larger degree in BA net- 
works. So hereinafter, all the judgements and operations 
are based on the degree of nodes. 

In order to enhance the network synchronizability, a 
few nodes with highest degree will be divided. Rank 
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FIG. 2: Behavior of the ratio of the eigenratio of network 
G(Np, m) to that of network G versus the number of divided 
nodes. As the number increases, the ratio is shown to reduce, 
leading to better synchronization. The average is taken over 
50 different network realizations. 



each node of a give network G according to its degree, 
the node that has highest degree is arranged in the top of 
the queue. Then, the network G(p, m) can be obtained by 
the following N p steps. First, carry out m-division on the 
top node in G, leading to the network G(jj, m). secondly, 
calculate all nodes' degree in G(-^,m), and rank each 
node according to it. Then, get the network G(jj,m) by 
dividing the top node in G(-^,m). Repeat this process 
Np times, with N p nodes have been divided in total, one 
will reach the network G(p, m). Since the randomicity is 
involved in dividing process, G(p, m) is not unique. In 
this letter, we focus on the difference between G(p, m) 
and G. 



IV. SIMULATIONS 

To explore how the structural perturbations affect net- 
work synchronizability, we compare the eigenratio r be- 
fore and after dividing processes. The BA networks of 
size N = 2000 and average degree (k) = 12 are used for 
simulations. In Fig. 2, we report the ratio R = r'/r 
against the number of nodes that are divided, where r 
is the eigenratio of the original network and r' after the 
operation. Here we set m = 3. With the probability p 
or the number of divided nodes increasing, the ratio R 
is observed to decrease, indicating the enhancement of 
synchronizability. In Fig. 2, it can be seen that to divide 
a very few nodes will sharply enhance the network syn- 
chronizability. R decreases to 0.7 when only 5 nodes are 
divided, and will drop to half after 0.6% nodes (i.e. 12 
nodes) are divided. To better understand the under- 
lying mechanism of synchronization and the reason why 



FIG. 3: The average distance L' and maximal degree k' max in 
G(p,m) vs p. L and fe max denote the average distance and 
maximal degree in original network G. We plot the relative 

changes 4- and ^ using squares and circles, respectively. 
One could see clearly, the dividing processes reduce the max- 
imal degree while increase the average distance. All the data 
are obtained by the average over 20 independent runs. 



these structural perturbations will greatly enhance net- 
work synchronizability, we investigate the behaviors of 
two extensively studied quantities, the average distance 
L and maximal degree /c max . In BA networks, the node 
with maximal degree is most probably the very node hav- 
ing maximal betweenness. As illustrated in Fig. 3, L will 
increase with p, while /c max will decrease. This results 
provide some evidences how the two factors affect the 
synchronization of systems. The maximal degree(i.e. the 
maximal betweenness) may play the main role in deter- 
mining network synchronizability. It is worthwhile to em- 
phasize that from the simulation results, we can not say 
anything about how the average distance affects the net- 
work synchronizability. L varies slightly, and probably 
has nonsignificant influence compared with the change 
of /c max - This results suggest that to reduce the maxi- 
mal betweenness of networks is a practical and effective 
approach to enhance the network synchronizability. 



V. CONCLUSION AND DISCUSSION 

Motivated by the practical requirement and theoretical 
interest, numbers of scientists begin to study how to en- 
hance the network synchronizability, especially for scale- 
free networks |4ll |42|. These methods keep the network 
topology unchanged, while add some weight into the sys- 
tem, thus the coupling matrix is changed. These excellent 
works do not need any new nodes, new edges or rewire- 
ment, but highly enhance the network synchronizability. 
However, it may require more intelligent system, in which 
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each node at least should remember its particular link- 
ing strengths. In this paper, we propose an approach to 
enhance the network synchronizability. This approach 
does not require any intelligence of nodes, but the net- 
work structure will be slightly changed. In some real- 
life communication networks such as Internet, the long 
length edge may cost much more than node and short 
length edge [431 l44j |. so if all the nodes x\, X2, ■ • ■ , £ m -i 
are in xq's vicinal locations, our method is feasible. 

Some recent works about network traffic dynamics re- 
veal that the communication ability of network, called 
network throughput |l lL IT^ . is mainly determined by the 
maximal betweenness, thus to steer clear of those hub 
nodes may enhance the network throughput 0> 01- This 
is just the case of network synchronization. Some method 
that can enhance the network throughput will enhance 



the network synchronizability too fllL l4ll 143 . l45 | . There- 
fore, we guess there may exist some common features 
between network traffic and network synchronization, al- 
though they seem completely irrelevant. We believe our 
work will enlighten readers on this object, and is also 
relevance to traffic control on networks. 
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